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ABSTRACT 

We investigate the use of two kinds of staggered fermion operators, smeared and unsmeared. The 
smeared operators extend over a 4 4 hypercube, and tend to have smaller perturbative corrections 
than the corresponding unsmeared operators. We use these operators to calculate kaon weak matrix 
elements on quenched ensembles at /3 = 6.0, 6.2 and 6.4. Extrapolating to the continuum limit, 
we find B K (NDR, 2 GeV) = 0.62 ± 0.02(stat) ± 0.02(syst). The systematic error is dominated by 
the uncertainty in the matching between lattice and continuum operators due to the truncation of 
perturbation theory at one- loop. We do not include any estimate of the errors due to quenching or 
to the use of degenerate s and d quarks. For the AI = 3/2 electromagnetic penguin operators we 
find Bf ,2) = 0.62 ±0.03 ±0.06 and Bf ,2) = 0.77± 0.04 ± 0.04. We also use the ratio of unsmeared 
to smeared operators to make a partially non-perturbative estimate of the renormalization of the 
quark mass for staggered fermions. We find that tadpole improved perturbation theory works well 
if the coupling is chosen to be a^(q* = I /a). 
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1 Introduction 



An important goal of lattice QCD is to provide reliable calculations of electroweak matrix 
elements. The major sources of error in present calculations are the use of finite lattice 
spacing, the use of one-loop perturbation theory to match continuum and lattice operators, 
and the use of the quenched approximation |l| . In this paper we address the first two errors 
for calculations using staggered fermions. In particular, we test the efficacy of "smeared" 
operators [0. These extend over a 4 4 hypercube, and thus are larger than the usual ("un- 
smeared") operators which are confined to a 2 4 hypercube. Nevertheless, in many cases they 
are closer to the continuum operators in the sense that the one-loop matching coefficients 
are closer to unity. 

We apply these operators to the study of three quantities. The first is the kaon B- 
parameter, Bk, which we study using both smeared and unsmeared operators. The initial 
motivation for introducing smeared operators was the discovery of large discretization errors 
in the results for unsmeared operators 0, |J. At the time, it was unclear whether the 
discretization errors were proportional to the lattice spacing a or to a 2 (up to logarithmic 
corrections). The smeared operators were designed to reduce possible 0(a) errors — they 
match onto continuum operators with no errors of 0(a) at tree level. It was subsequently 
realized that the 0(a) parts of the lattice operators do not contribute to B K , and that the 
errors are automatically of 0(a 2 ) for both the unsmeared and smeared operators || |6|. [7|. 
This theoretical argument has since been tested numerically [0, |S|. Thus the interest in 
using the smeared operators is that they provide an estimate of the error in the matching of 
continuum and lattice operators. The results from different operators should differ at finite 
a, but agree upon extrapolation to a = 0, up to higher order perturbative corrections. 

Preliminary results from this study of Bk were presented in Ref. ||, and one of our 
purposes here is to present final results. Although, by present standards, these come from 
a small statistical sample, the errors are nevertheless small enough to assess the impact of 
smeared operators. We have also improved our estimates of the error due to the truncation 
of the perturbative matching factors, using the method introduced in Ref. 0. 

Our second application is the calculation of B^ 2 and B^ 2 . These B-parameters, and 
in particular B^ 2 , determine the size of the electromagnetic penguin contribution to e'. In 
contrast to Bk, the use of staggered fermions for these quantities offers no clear advantage 
over Wilson fermions. However, since the systematic errors in the results with the two types 
of fermion are different, an important check of the reliability of the lattice calculations is to 
show that the two formulations give consistent results in the continuum limit. To this end 
we compare the staggered data with the recent results obtained using Wilson fermions at 
(3 = 6/g 2 = 6 in the quenched approximation ||. 

The calculation of B% and B^ 2 demonstrates the importance of using several discretiza- 
tions of continuum operators. It turns out that one cannot use the unsmeared operators 
because the one-loop correction to the matching coefficients approaches 100% ]TU| . On the 
other hand, the one-loop corrections are much smaller for the smeared operators (~ 25%), 
and we can use them to calculate Bj and B^ 2 . The only way in which one could use 
the unsmeared operators would be to develop a non-perturbative method of determining the 



matching coefficients (which is possible in principle using external quark states jjTT 



Our final application concerns the calculation of the matching relation between the lat- 
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tice and continuum regularization schemes, particularly in cases where the reliability of 
perturbative estimates is questionable. The ratio of such matching factors for two differ- 
ent discretizations of an operator can be estimated non-perturbatively by taking ratios of 
appropriate matrix elements. The non-perturbative results so obtained can be used to test 
the reliability of the one-loop perturbative estimates. In particular, one can use the results 
to fix the scale q* at which to evaluate the coupling constant entering into the perturbative 
expressions. For the pseudoscalar density we find that q* ~ 1/a. Our conclusions are, how- 
ever, preliminary, since we do not have results at enough values of lattice spacing to check 
the extrapolations we use to remove discretization errors. 

We use this method to assess the reliability of the matching factor, Z m . Z m relates the 
bare lattice mass to the continuum mass in, say, the MS scheme, and is a crucial ingredient 
in the calculation of continuum light quark masses from the lattice. Recent work has sug- 
gested that continuum quark masses are smaller than previously thought, but this is based 
on trusting one- loop perturbation theory for Z m |12| . For staggered fermions, the one-loop 
contribution to Z m is large, roughly a 60% correction at (3 = 6, even after tadpole improve- 
ment [0]. This casts doubt on the reliability of the perturbative Z m , and therefore also on 
the extracted value of quark masses. We find, however, that our partly non-perturbative 
estimate suggests that the one-loop evaluation is close to the correct answer. 

The organization of this paper is as follows. In the following section we describe the 
method we use to match continuum operators to lattice operators composed of staggered 
fermions. In sec. B] we give a short description of the numerical methods and data sample. 

3/2 

The three subsequent sections contain our analysis and results for Bk, B 7 8 , and the non- 
perturbative ratios of matching factors, respectively. We close with some conclusions. 



2 Theoretical Review 



In this section we explain our method of calculating 5-parameters using staggered fermions. 
This requires combining a variety of results already in the literature, and we focus here only 



on the essential details. For a more extensive description of the method see Ref. |13| 
The continuum operators of interest are 



Qk 

,3/2 



'A' 

Sal^db [Ubl^U a - dbrfda ] + s a ^u b u b j^Rd a • 



(1) 

(2) 
(3) 



where a, b are color indices and 7^' = 7^(1 ± 75). All operators are in Euclidean space, 
and we use hermitian Gamma matrices. The superscripts on Q 7 8 indicate that these are the 
I = 3/2 parts of the operators ^7,8, i-e. the I = 1/2 component has been removed. We make 
this restriction because the calculation of the matrix elements of the I = 3/2 parts is much 
simpler. In the limit of exact flavor SU(3), which is the limit we work in here, the / = 3/2 
parts give rise only to "eight" diagrams, i.e. those in which the quark fields in the operator 
are contracted with fields in the external mesons. These are the same type of diagrams which 
contribute to the matrix element of Qk- The 7=1/2 parts, by contrast, give rise also to 
"penguin" or "eye" diagrams, which are much more difficult to calculate. The restriction to 
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the I = 3/2 parts of Qr^ does not, however, diminish the phenomenological interest in the 
results, because 8 2 are the only operators which give an imaginary part to the K + — > 7r + 7r° 
amplitude. 

To form P-parameters we also need the matrix elements of the axial and pseudoscalar 
densities, 

A^ = s 7 M 7 5 d , P = Sa^da . (4) 

We can then define 





] \Qk\ 


K°) 






(8/3)(K°| 


A 4 \0) (0\A 4 \K°) 
(tt + \Q 3 7 /2 \ 


(2/3)(^°| 


P\0) (0\P\K°) - 


(K°|A 4 |0) (0|A 4 |K°) 



8 ~ 2(iT°|P|0) (0|P|K°) - (1/3)(AT |A 4 |0) (0|A 4 |i?°) ' 1 ; 

All external particles have been assumed to be at rest. For brevity, we have used SU (3) flavor 
symmetry to rewrite all the denominators in terms of kaon matrix elements. Note that, in 
general, both numerators and denominators of these ratios depend upon the renormalization 
scale fj, and the scheme used to define the operators. When quoting physical values we use 
the NDR scheme, i.e. MS renormalization combined with a particular set of rules for treating 
7 5 away from four dimensions, and choose the renormalization scale to be \i = 2 GeV. 

The extraction of the above matrix elements using staggered fermions is complicated by 
the mixing between the spin and flavor degrees of freedom. As explained in Refs. |L0], [13| , 
we proceed in two stages. We first match the continuum matrix elements onto those in 
an "enlarged" continuum theory, and then match from that theory onto the lattice. The 
enlarged theory differs from QCD by having eight copies of each physical quark. The eight 
copies of the strange quark spinor are collected into two 4x4 matrices Spj, an d S'a b , where f3 
is a spinor index, and b = 1 — 4 is a "staggered-flavor" index. Similar fields are constructed 
for the up and down quarks. The correspondence between matrix elements in the continuum 
and in the enlarged theory is 

(0\A 4 \K°) = yX(0|£( 74 75 ® &D\K° G ) , (8) 

(K°\A 4 \0) = ^(K'£\S'(w 5 ® £ 5 )£>'|0> , (9) 

(0\P\K°) = </^-<0|S(7 5 ® &D\K° G ) , (10) 



N 



f 



(k°\p\o) = ^(K'Sis'i^ ® £ 5 )r/|o> , (ii) 

(K°\Q K \K°) = 1-(K%\ [(V - A) x P\j + [(V - A) x P\ n \K%) , (12) 



N f 
1 



(n + \Q 3 7 /2 \K + ) = ±(K'g\ 2 [(P - S) x P], + [(V + A) x P] u \K G ) , (13) 
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(TT+IQf \K+) = -L(Kg\ 2[(P - 5) x P] 7/ + [(V + A) x P] 7 \K%) . (14) 

We have used flavor symmetry to rewrite all matrix elements in terms of those between 
external kaon states. The notation for matrix elements in the enlarged theory is that of 
Refs. m |10| , and we give only a brief summary. The matrices appearing in bilinears are 
tensor products of spin and staggered-flavor matrices. For example, (75 ® £5) indicates 
a pseudoscalar density with staggered-flavor matrix 75. The states \Kq) and \Kq) are 
those created from the vacuum by D(j 5 x £ 5 )S and D'(j 5 x ^)S', respectively. The states 
are normalized, which leads to the factors of Nf = 4, the usual multiplicity factor for 
staggered fermions. The subscript G indicates that these states are the pseudo-Goldstone 
bosons corresponding to the axial U(l) symmetry which is unbroken when one discretizes 
the enlarged theory using staggered fermions and takes the chiral limit. The notation for 



four-fermion operators is from Ref. |T0[. For example, [(V — A) x P\j represents the one 
color loop contraction of the four-fermion operator with spin structure 7^ ■ 7 M — 7^75 ■ 757^, 
and in which both bilinears have staggered flavor 75. Finally, the factor of 2 on the r.h.s. 
of Eq. (|I2|) arises from the difference in the number of Wick contractions in QCD and the 
enlarged theory due to the use of "primed" quarks in the latter theory. 

All these equalities hold identically between the quenched versions of the two theories. 
In the presence of internal fermion loops, they hold if, in the enlarged theory, each loop is 
multiplied by 1/8, i.e. the fermion determinant is taken to the power 1/8. 

The next step is to relate the operators in the enlarged continuum theory to those in the 
corresponding lattice theory. The enlarged theory has been chosen so that it is the continuum 
limit of a discretized theory in which one uses a single staggered species for each 4x4 matrix 
field, e.g. xs for S and xs' for 5". This means that the relation between operators is simple 
at tree level. For example, using the notation of Ref. 0, 



S(75 ® &>)£> = Xs(7s ® &Xd [1 + 0(a) + 0(a s )} . (15) 



Here (75 (g) £5) is that matrix in the space of possible positions in a 2 4 hypercube which 
corresponds to the continuum spin-flavor matrix (75 ® £5). Similar relations hold for the four 



fermion operators |TU 



In this paper we use one-loop perturbation theory to match operators in the enlarged 
continuum theory to those on the lattice. In other words, we include the terms of 0(a s ) in 
equations such as ([15]). More precisely, we use the "horizontal matching" procedure discussed 
in detail in Ref. ||. This consists of two steps. We first match the lattice and continuum 
operators at an intermediate scale q* ~ 1/a using 

)cont(„*\ _ f^lat 1 a Usi a ) \ " ( i „..\f1 lat 

Ait 



0?* if) = OT + -^f^ £ (- 7 §V Mq*a) + c %3 ) Of + 0(a 2 ) + 0(a) . (16) 



Here O 1 ^ 1 are bare lattice operators, 7*- - 1 is the one-loop anomalous dimension matrix, and 
the Cjj are the one-loop matching coefficients. We give, below, numerical values for the c^- of 
interest . The second step is to evolve the result from q* to the final scale \x = 2 Ge V using 
the continuum two-loop anomalous dimension matrix. The two-loop anomalous dimensions 
for the continuum operators of interest are collected in Ref. 0, and we do not repeat them 
here. 
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The point of this procedure is to account for the fact that there are two scales in the 
problem: the matching scale q* ~ 1/a, and the final scale \i. These can differ substantially — 
indeed the ratio q*/fi can be as large as 10 in our calculation. The coupling constant 
can thus be quite different at the two scales, and it is important to make sure that the 
appropriate coupling is used at each stage. Furthermore, it becomes necessary to sum up 
the leading logarithms of q* / //. Both of these requirements are accomplished by horizontal 
matching. The scale in the coupling in the first step should be q*, and the renormalization 
group evolution in the second step sums up the leading logarithms. Of course, since we 
truncate perturbation theory, there are errors in both steps. In the matching relation ([16]) 
the truncation errors are of 0[a s (q*) 2 ], while the errors in the continuum evolution are 
proportional to both 0[a s (q*) 2 } and 0[a s (fi) 2 ]. We note, however, that the errors from the 
continuum evolution are the same for any choice of lattice discretization of the continuum 
operator. This point will be important in our discussion of results for B^ using different 
lattice operators.^ 

To use horizontal matching we need to choose the intermediate matching scale q*. If one 
worked to all orders in perturbation theory, the final result would be independent of q*. This 
is not true if one truncates perturbation theory: different choices lead to results differing by 
0[a(q*) 2 ]. One should choose q* to be, roughly speaking, the average momentum flowing 



through the quark-gluon vertices in the matching calculation [16|]. It is difficult, however, 
to make this into a precise prescription for matching calculations involving operators with 
non-zero anomalous dimensions. The only method we know of was proposed in Ref. 0, and 
involves applying the BLM prescription of Ref. [T7| to the calculation of the two-loop lattice 
anomalous dimension matrix. This calculation has not been done, and so we have chosen 
to use a range of values for q*. Since the lattice and continuum operators are constructed 
to have the same matrix elements at long distances, the dominant contributions to the 
matching calculation are from momenta near the lattice cut-off (as long as the continuum 
renormalization point is also of this size). Thus we take q* = K/a with K a constant. Since 
we use tadpole improved operators we expect that K « 1 rather than K m it, as explained 
in Ref. ||16|| . So we have chosen q* = 1/a for our central values, and used q* = it /a in order 
to estimate the uncertainty due to the truncation of perturbation theory^ This is certainly 
a crude estimate, but it is the best that we can do. 

In our preliminary analysis of Bk we did not use horizontal matching, but rather used 
the one-loop matching relation Eq. ( |T6"D connect directly from the lattice to the final scale 
jj, ||. To estimate the truncation error we took the difference between the results of using 
a s (n) and a s (q*) in the matching equation. Our present methods both of matching and of 
estimating the truncation error are more reliable, although, as we will see below, the final 
answer is little changed. 

The MS coupling constant appearing in Eq. ( [IB] ) is determined using the method of Ref. 

4 For further discussion of horizontal matching, and in particular of its relation to the exact matching 
formula of Ji ]T^|, see Ref. || . 

5 Reference |16| advocates the use of a different coupling constant, ay, rather than a-^g. ay is defined in 
terms of the quark-antiquark potential. We prefer to use a-^r because this is the scheme used in continuum 
calculations of coefficient functions, and the matching formulae are simpler if one uses the same coupling in 
the continuum and on the lattice. The two couplings are quite similar — our choice of q* = 1/a in the MS 
scheme corresponds to using ay(1.6/a). This is indeed a typical value for q* for tadpole improved quantities. 
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16], which incorporates tadpole improvement. We first solve 



47T 

- log □ = —a v {3Al/a) [1 - 1.185ay(3.41/a)] (17) 
o 

for ay, where □ is the expectation value of the plaquette normalized to be unity in the 
continuum limit. We then convert to the MS scheme using 

a m (3Al/a) = a v (e 5/G 3.41/a)(l + 2a v /n) . (18) 

This coupling is then run to other scales using the two-loop /3-function. 

The lattice operators we use are of two types: the standard operators, which are contained 
in 2 4 hypercubes (and which we refer to as "unsmeared" operators); and the "smeared" 
operators introduced in Ref. B, which live on 4 4 hypercubes. The smeared operators have 
the same form as the unsmeared, except for the replacement 

x(y + A) - W x(y + a + 2A[i - , (19) 

where y is the position of the origin of the 2 4 hypercube, and A is the position within the 
hypercube. This definition makes the average position of the quark field lie at the center of 
the hypercube, which is why there are no 0(a) corrections when matching to the continuum 
operator at tree level. Both types of operator are made gauge invariant by fixing to lattice 
Landau gauge. Both are also tadpole improved, using the mean link determined from the 
plaquette, Uq — □. 

Although the discretization errors in the matching equation are, in general, of 0(a), it 
turns out that for the matrix elements of interest the corrections are actually of 0(a 2 ). This 
is because the 0(a) parts have the wrong staggered-flavor |5|, |TJJ, and by using staggered 



flavor £ 5 for the incoming and outgoing states we project against them up to errors of 0(a 2 ). 
For the same reason, we need only include in the matching equations those operators which 
have the correct staggered-flavor. This is fortunate, since typically many operators of the 
wrong flavor are needed to match the continuum operators. 

The matching matrices Cjj for all continuum operators of interest, and for both smeared 



and unsmeared lattice operators, can be determined from the results of Ref. ||10|| . When 



doing this one must be careful to account for the following three issues. First, the numerical 
results for the q,- are quoted in Ref. fLOl for continuum operators defined in the DREZ , 



rather than the NDR, scheme. Thus one must convert between these schemes, and a recipe 
for doing so is given in [|K|. Second, tadpole improvement in [T(J is based on the average 



trace of the link in Landau gauge, rather than on the average plaquette which we use here. 
This causes a small change in the coefficients 

c i:j = Cij(Ref. [10]) + iV B (4/3)(7r 2 - 9.17479)^ , (20) 

where Nb = 1 for bilinears and Nb = 2 for four fermion operators. Third, the coefficient 
of the logarithm in the matching equation (|T6|) has been changed from q*a/ir to q*a. This 
shifts the Cij of Ref. [|K| by 7^ In tt. This change allows one to see the size of the matching 
corrections for our standard value q*a = 1 directly from the size of the c^-. 
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As the procedure for calculating the matching coefficients is rather involved, we collect 
here the relevant parts of the final results. For the bilinears, the matching involves no mixing, 
only multiplicative renormalization. The results are 



C A 

UN 
Cp 



0.9264, 
-39.1414 



C A 

„SM _ 
Cp - 



= -1.0120. 
-8.8882. 



(21) 
(22) 



Here ca is a shorthand for the diagonal coefficient for the axial density, etc., and the super- 
script refers to whether or not the quark fields are smeared. These coefficients are multiplied 
by aj^g(g*)/47r, which, in our calculation, varies in the range 0.01 — 0.015. Thus these cor- 
rections are numerically small for the axial currents, (~ 1%), and of moderate size for the 
smeared pseudoscalar density (~ 10%). For the unsmeared pseudoscalar density, however, 
they are large (~ 50%), suggesting that two-loop terms will be important. This is an example 
of the advantage of using smeared operators. 

For the four-fermion operators of interest here, we can decompose the matrices of match- 
ing coefficients into two parts. The first is a square matrix with indices running over operators 



i, j = ( [(V - A) x P] Jt [(V - A) x P] Ih [{V + A) x P) h [{V + A) x P] H ) 
for which we find 



(23) 



SM 



^-4.3079 0.4611 0.2638 -0.7913^ 
-1.0002 0.0761 -2.1102 
0.2638 -0.7913 -8.9740 -1.5387 
-2.1102 -4.9998 1.4093/ 

-2.1376 -5.4068 0.0850 -0.2551^ 
-5.8701 -0.7474 -0.6802 
0.0850 -0.2551 4.5595 -2.5366 
V -0.6802 -0.1299 -2.6608/ 



(24) 



(25) 



These coefficients are needed for all three 5-parameters, and the corrections are of moderate 
size for both smeared and unsmeared operators. The second matrix we need is rectangular, 
having indices^ 



i 

J 



[{P 
i(P 



S)xP]j, [(P-S)xP] n ) , 

S) x P] h [(P -S)x P] n , [(P + S)x P]j, [(P + S)x P] n ) 



For this we find 



0.7562 
-3 

-2.8862 
-3 

,3/2 



-14.7309 
-34.4364 

-3.1611 
-3.3692 

?3/2 



5.3559 -16.0674 \ 
-42.8464 j 

1.6759 -5.0277 \ 
-13.4072 ) ' 



(26) 
(27) 



(28) 
(29) 



These results are needed for B 7 and B 8 , and it is here that the use of smeared operators 
is most important. In particular, the numbers in the second row in Eq. fl29|) are much smaller 



3 Note that, at one-loop, tensor operators do not appear. 
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(3 6.0 6.2 6.4 

Lattice Size 24 3 x 40 32 3 x 48 32 3 x 48 

Number of lattices 13 23 24 

Samples per lattice 12 2 

Quark masses 0.01, 0.02, 0.03 0.005, 0.01, 0.015 0.005, 0.01, 0.015 

l/a in GeV from m p 1.9 2.6 3.45 

Minimum M n in GeV 0.454 0.376 0.414 

Average plaquette 0.5937 0.6137 0.6306 

a m (n/a) 0.134 0.125 0.117 

a m (l/a) 0.192 0.173 0.158 

^(2 GeV) 0.188 0.190 0.191 

l/a for a m (2 GeV) = 0.190 1.94 2.6 3.4 



Table 1: Parameters of numerical simulations used in this analysis. 



than in the corresponding row in Eq. (^). Much of this reduction can be traced back to the 
similarly large difference in the coefficients cp shown in Eq. (p2|). For the range of couplings 
we study, it turns out that the one-loop matching correction for the operators appearing in 
the numerator of B^ 2 and B\' 2 actually exceeds 100% if one uses the unsmeared operators 
with q* fa l/a. Our results for these quantities are, therefore, obtained exclusively with 
smeared operators. 



3 Numerical details 

Our results are based on ensembles of lattices at three different lattice spacings. A summary 
of the important parameters is given in Table [I]. The lattices at (3 = 6 are part of a 
sample previously used to study the hadron spectrum, and a complete description of how 
we generated the lattices and calculated quark propagators can be found in Ref. flI8| . The 
lattices at (3 = 6.2 and 6.4 have been discussed previously in Refs. || |19| . They were 
generated using overrelaxed and Metropolis sweeps in a 4:1 ratio, and separated by 1000 
sweeps. They were divided roughly equally into two independent streams at (3 = 6.2, and 
three at j3 = 6.4. Quark propagators were calculated using the conjugate gradient algorithm. 

We consider here only pions composed of two degenerate quarks. For the quark masses 
that we use, the masses of the lattice pions bracket that of the physical kaon at each (3. To 
illustrate this we include in the table the mass of the lightest pion, converted to physical units 
using l/a determined from m p . These scales, also listed in the table, have been previously 



reported by one of us from a fit to the hadron spectrum [19]. We refer to these below as the 



'm p scales'' 



We also give the values of needed in the one-loop matching between lattice and 
continuum operators. These are obtained using Eqs. ([H]) and ( |18D from the plaquette values 
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listed in the table. Since in the end we run our results to 2 GeV, a consistency check on 
the applicability of tadpole improved perturbation theory is that the values of 0^(2 GeV) 
agree for different (3. The table shows that they are in reasonable agreement, although there 
is a small systematic increase as a decreases. 

This systematic dependence on a suggests using a different set of scales, defined so that 
a^g(2 GeV) is independent of a. We choose «ms(2 GeV) = 0.190, and the resulting scales 
are listed in the last row of the table. We refer to these below as the "a scales". We use 
these scales to determine our central values for matrix elements, and use the m p scales to 
estimate the error introduced by the uncertainty in determining a. 

Our method for calculating 5-parameters was developed for calculating B^ with stag- 
gered fermions p0| . It has also been used for Wilson fermions pi . A detailed description 



of the latter application, and the extension to other operators is given in Ref. ||. Here we 
confine ourselves to a brief description of the features which are special to the present work. 

The essential feature of the method is the use of wall sources (i.e. sources confined to 
a single time-slice and extended over the entire timeslice), which are designed to create 
particles with specific quantum numbers. In the current applications we want to create 
only the lattice pseudo-Goldstone pion at rest, and this can be accomplished using a linear 
combination of two wall sources, as explained in Ref. |18|| . The only other states created are 
rho mesons and excited pions (tt'). Both of these are, however, considerably more massive 
than the pseudo-Goldstone pion, and their contribution can be largely removed by moving 
far enough away from the source in Euclidean time. 

The basic method is then to calculate ratios such as 

K{> S(W(t 1 )Ey'A,(y',t))(E f/ A,(y",t)W(t 2 ))- { ] 

Here the operators are the lattice versions of the continuum operators obtained after the 
matching explained in the previous section has been carried out. The wall sources are denoted 
by W and are located at times t\ and t 2 . The operators are placed at an intermediate time 
t satisfying ti <^ t <^ t 2 , and are summed over the V y spatial hypercubes. Finally, the 
expectation values in Eq. ([30|) are averages over quenched configurations. 

The expression in Eq. ( p0[ ) is designed so that for large enough t — t\ and t% — t it is 
independent of t, and gives Bk directly. The point is that the exponential factors from 
Euclidean time evolution cancel if only a single state, here the pseudo-Goldstone pion, con- 
tributes. This has to be true separately for the numerator and denominator, i.e. both must 
exhibit a "plateau" over intermediate times in which they are independent of t. Thus the 
signal can be improved by averaging the numerator and denominator over the time-slices 
in this plateau. This is what we do in practice, using the same range of times for both 
numerator and denominator. 

This method has several positive features. First, it involves no fitting. Second, statistical 
errors are reduced by directly calculating 5-parameters rather than the matrix elements 
themselves. Third, the ability to average over all the spatial points and over a number of 
time slices improves the statistics. This is possible because we know the propagator from the 
wall sources to all points in the lattice. For the same reason, we can calculate matrix elements 
of non-local operators without additional quark propagators. And, finally, we can use the 
same set of quark propagators to calculate matrix elements with different choices of lattice 



9 



operator. The latter two features are particularly important for staggered fermions where 
most operators are non-local. If we placed the source of the propagators at the position of 
the operators we would need 2 4 sources for unsmeared operators (4 4 for smeared operators). 

To calculate statistical errors, we use single elimination jackknife in the following manner. 
On each jackknife sample, we first match at the scale q*, then linearly interpolate the results 
to physical kaon mass, and then evolve to the final scale /i = 2 GeV, and, finally, take the 
ratio which defines the 5-parameter. The error is then obtained from the variation between 
samples. 

We end this section with some details specific to our particular lattices. 

• At (5 = 6, we use wall sources in Coulomb gauge, while on the other lattices the sources 
are in Landau gauge. The choice of gauge should have no impact on the final result 
(as long as the non-locality introduced by the gauge fixing does not extend from the 
source to the operator JI8| ) but might affect the statistical errors. 

• At p = 6, we use periodic boundary conditions (PBC) in space, and Dirichlet bounday 
conditions in time. We place the wall sources next to the boundary, i.e. at t% = and 
t 2 = 39. The merits of this procedure have been discussed in Ref. We only note 
here that we use a plateau region of t — 10 — 29. 

• The lattices at /3 — 6.2 and 6.4 have too short an extent in time to follow the method 
adopted at /3 = 6. Instead we use PBC in all four directions, having first periodically 
doubled the lattice in time. We place wall sources at t\ — and t 2 = 72, and use the 
plateau region t = 32 — 40. Note that the propagator with source at t\ — it 2 = 72) 
is, due to the PBC and doubling of the lattice, the same as if the source was at t\ = 48 
(i 2 = 24). Because of this, and the fact that each source produces pions propagating in 
both forward and backward directions, we can make a second measurement with the 
plateau region at times t = 8 — 16. We treat these two results as independent in our 
jackknife error analysis. 

To study possible sources of systematic errors due to contamination from excited states 
or "off-shell" matrix elements we have also made the following measurements. Taking 
the same two wall sources to lie at t% — and t 2 = 24 we use the region t = 32 — 40 
to obtain an estimate of the "off-shell" matrix elements, i.e. those in which both of 
the pions approach the operator from the same side. Placing the sources at ti — and 
t 2 = 48, and considering the plateau region at t — 23 — 24, we get an estimate of the on- 
shell matrix elements having a greater contamination from excited states, but a smaller 
contribution from off-shell matrix elements. Further discussion of these constructions 
and the associated systematic errors is given in Ref. ||. In practice, we find that the 
various methods yield very similar results for the ^-parameters we consider, and that 
these sources of error are considerably smaller than others as discussed later.Q Thus 
we give no further details. 



7 The same is not true for the auxiliary parameters Ba and By defined in Ref. 0] . For these we have no 
useful results at ft = 6.2 and 6.4. 
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Operator 


Q* 


(3 = 6.0 


(3 = 6.2 


(3 = 6.4 


a = 


UN 


1/a 

I 


0.728(9) 


0.677(9) 


0.685(24) 


0.628(20) 


UN 


n/a 


0.725(9) 


0.669(9) 


0.677(24) 


0.614(20) 


SM 


1/a 


0.681(6) 


0.647(6) 




0.604(15) 


SM 


71/a 


0.694(6) 


0.655(6) 




0.606(15) 


Ratio UN/SM 


1/a 


1.069(11) 


1.046(13) 




1.018(32) 


Ratio UN/SM 


Ti j a 


1.044(11) 


1.021(14) 




0.991(34) 



Table 2: Results for S^(NDR, 2 GeV) at the physical kaon mass, using the a scales. "UN" 
and "SM" refer to unsmeared and smeared operators, respectively. 



4 Results for Bk 

Using the methods described in the previous two sections we extract the continuum Bk at 
li — 2 GeV for each lattice spacing. The only feature of the analysis not discussed above is 
the interpolation to the physical kaon mass. This we do by fitting Bk itself (for our three 
mass points) to a linear function of the squared lattice kaon mass, m 2 K lat . This is reasonable 
for staggered fermions because chiral symmetry constrains Bk to have the same form as 
in the continuum and in particular to be finite in the chiral limit [O, E3. For degenerate 
quarks, the explicit form is 



B K = B 



l-3y\ny + by + 0(y 2 ) , (31) 



where y = m^-/(47r/ 7r ) 2 is the usual chiral expansion parameter. For the our range of m Klat 
the y\ny contribution is well represented by a linear function. We note in passing that 
for Wilson fermions there is, in general, an additional term proportional to a/y in Eq. ( |31"D 
because of the explicit breaking of chiral symmetry. 

We present results only for the a scales, i.e. the lattice spacings determined by requiring 
that 0^(2 GeV) = 0.190. These results are collected in Table 0, and displayed in Fig. |I|. We 
include the results of an extrapolation to the continuum limit assuming quadratic dependence 
on a. Note that we do not have results for smeared operators at (3 = 6.4. Our results are not 
extensive enough either to test whether the dependence on a is indeed quadratic, or whether 
terms of higher order than quadratic are needed for our range of lattice spacings. The best 
confirmation of the validity of the quadratic dependence comes from the work of the JLQCD 
collaboration, who have more extensive results than ours for both the unsmeared operator 
and its gauge-invariant version JF|. 

We can use our results to estimate the systematic errors arising from the choice of lattice 
spacings, the truncation of perturbation theory in the matching factors, and the contamina- 
tion from excited states and off-shell matrix elements. 

We begin with the dependence on the choice of lattice spacings. This we estimate by 
comparing the results using the a scales (listed in Table ||) to those obtained using the scales 
determined from m p . The values of Bk at each lattice spacing change by no more than 
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0.002, and the maximum change in the extrapolated value is 0.004. Thus we take ±0.004 
for our estimate of this systematic error. This is much smaller than the statistical errors. 

The error due to the truncation of perturbation theory in the matching factors can be 
estimated in two ways. The first, discussed in sec. 0, uses the q* dependence of the result. 
We take the error to be the difference between the extrapolated results with q*a = 1 and 
q*a = 7i. The table shows that this is comparable to the statistical errors. The largest 
difference is that for the unsmeared operators, and we take this for our estimate, yielding 
±0.014. This 3% error is a reasonable estimate for a two-loop correction given that the 
one-loop matching corrections for unsmeared and smeared operators are ~ 5 — 15%. Our 
estimate of this error turns out to be the same as that quoted in our preliminary result 0, 
although the method we use here is more reliable, as explained in sec. ||. 

The second way of estimating the perturbative error is to compare the results using the 
two types of operator. If the matching factors were correct then they should yield the same 
result in the continuum limit. We estimate the error as half the difference between the 
smeared and unsmeared results at q*a = 1, and thus obtain ±0.012. In fact, this is likely to 
be an overestimate of the difference in the results from the two operators. This is because 
there is an extra data point at (3 = 6.4 for the unsmeared operators, and, as can be seen from 
Fig. 0, this shifts the extrapolated result away from that of the smeared operator. Perhaps a 
better way of estimating the error is to take the ratio of the results using the two operators, 
and extrapolate this to the continuum limit. This removes the data point at /3 — 6.4 (since 
we have no smeared result there), and also accounts for the correlations between the results 
for the two operators. The results for this ratio are given in Table and show that although 
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the results from the two operators are significantly different at finite a, they are consistent 
in the continuum limit. Because of this, we take our estimate of the perturbative error from 
the dependence on q*. 

Finally, we discuss the errors due to contamination in the matrix elements from excited 
states and off-shell contributions. The only excited state which is allowed to contribute by 
the symmetries is the tt', i.e. the radially excited pion. There are no contributions from p 
mesons (a point not realized in Ref. ||). As mentioned in Section the size of the first 
effect can be estimated by comparing the results with the two sets of sources (ti = 0, t 2 = 72 
versus t\ = 0, t 2 = 48), while off-shell contamination can be estimated by studying B K {t) in 
regions where the off-shell contributions dominate. We find that the combined shift in Bk, 
after averaging the results from the two operators, and extrapolating to the continuum limit, 
is ~ —0.003. Since this estimate is approximate and small, we do not include this shift in 
our final result, but instead include it as part of the above overall systematic error. 

Putting this all together, we can now quote our final result. For the central value we 
use q* = 1/a, and take the average of the results from the two operators. We use the larger 
of the statistical errors (that for the unsmeared operators). And we estimate the overall 
systematic error by combining linearly those from the choice of lattice scales (0.004), from 
the truncation of perturbation theory (0.014), and from the contaminations (0.003). Thus 
we quote 

B K (NDR, 2 GeV) = 0.62 ± 0.02(stat) ± 0.02(syst) . (32) 

This result is consistent with our preliminary number quoted in Ref. |5j (0.616±0.020±0.017), 
although the precise agreement is somewhat fortuitous given that our method of matching 
has been improved. 

Our results for the unsmeared operator can be checked by comparing them to those from 
the JLQCD collaboration 0. They have results at /3 — 6.0, 6.2, and 6.4, on lattices of 
the same spatial sizes as ours, but with statistical errors two or three times smaller. They 
also have results on larger lattices at /3 — 6 and 6.4, and at smaller values of (3. A direct 
comparison is possible because they use a method of matching and determining a-^ which is 
very close to ours if we set q* = 1/a. At /3 = 6 and 6.2, our results are larger by about 0.025, a 
two standard deviation difference. At /3 = 6.4 our number is consistent with theirs from a 32 3 
spatial lattice. They find, however, that Bk decreases on larger lattices, suggesting that our 
number at /3 = 6.4 may be afflicted by finite size errors. Nevertheless, our extrapolated value 
is only 1.5 standard deviations (i.e. 0.03) above theirs. This comparison gives us confidence 
that our results are correct within the quoted errors, and in particular that our procedure of 
doubling the lattice in the time direction has not introduced additional systematic errors. 

The JLQCD collaboration has also used the gauge-invariant version of the unsmeared 
operator. Averaging this with the Landau gauge operator, they quote a preliminary result 
5^(NDR, 2 GeV) = 0.587 ± 0.007(stat) ±0.017(syst), or adding errors in quadrature: B K = 
0.59 ±.02. The agreement with our result of Bk = 0.62 ±.03 is gratifying, since the operators 
we use have entirely different perturbative and power corrections.^ 

We note that the both calculations are systematics limited, and that of the systematic 
errors, the most important one quoted by JLQCD is estimated by comparing results for 

8 We can improve the agreement by dropping our unsmeared data point at = 6.4, which, as noted above, 
might be afflicted with finite size errors. 
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different operators. In their preliminary report, JLQCD ascribes such differences to terms 
of order a^g(/i) 2 , where \i = 2 GeV is the conventional scale at which the answer is quoted. 
On this point we disagree in principle. We certainly agree that errors of order a^(^) n 
are introduced when one uses n-loop evolution to the final scale, but these corrections are 
universal and should not appear as differences between lattice operators. If, for example, we 
compare our one-loop corrected smeared and unsmeared operators, the connection between 
the two is 

QT(») = QTM {l + 0[a m (q*f] + 0(a 2 )} , (33) 

where we have taken the same q* = K/a in renormalizing both operators.0 Since aj^r(K/a) 2 
vanishes as a — > 0, we conclude that when correctly extrapolated, the operators should give 
the same result at a = 0. Of course one would need rather precise data to make a fit including 
the ay^(K/a) 2 term, but in principle it could be done if more precision were required. If 
one uses a simple a 2 extrapolation, however, the a 2 term in Eq. (^) will not extrapolate to 
zero, but instead to an artifact of size «ms(?*) 2 - 



5 Results for B 7 and B 8 

We have evaluated B^ 2 and B^ 2 using almost the same method as for Br. The only 
difference concerns the extrapolation to the physical kaon mass. This we have done separately 
for matrix elements appearing in the numerator and denominator of the definitions Eqs. (^) 
and ©, prior to the evolution from q* to \i = 2 GeV. 

Our results are summarized in Table ^. The most striking feature is the very strong q* 
dependence of the results for unsmeared operators. Indeed, as one goes from q* = n/a to 1/a, 
which causes ajj^i^*) ^° increase by roughly 40%, Bj ' 8 2 change sign because the negative one- 
loop matching contribution exceeds the tree-level contribution. Clearly, we cannot use one- 
loop matching for the unsmeared operators. We stress that the large perturbative corrections 
for unsmeared operators do not invalidate the results for smeared operators. There will 
always be choices of discretization procedure for which the perturbative series is poorly 
convergent at the lattice spacing one is working. 

For completeness we mention that we expect a similar problem to render the gauge- 
invariant unsmeared operators unsuitable for a calculation of B^% . The dominant contribu- 
tion to the matrix elements comes from the PP part of the operator. Since this part of the 
operator is local, and does not require gauge links, it is the same as the corresponding part 
of the unsmeared operator. But it is this part of the unsmeared operator which leads to 
the bulk of the large one-loop matching corrections seen in sec. [| Although in principle it 
is possible these large corrections will be canceled by the as yet uncalculated contributions 
from the other parts of the gauge-invariant operator, this seems unlikely in practice. 

The situation is much improved for the smeared operators — the q* dependence, while 
more significant than for Bk, is only at the 10% level. This is a much larger uncertainty 
than the statistical errors or that due to the choice of lattice spacings, or that from the 

9 Using a different choice of q* for the two operators would lead to an additional factor coming from the 
evolution between the two scales. However, since this factor tends to unity in the continuum limit, our 
conclusion is unaltered. 
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Operator 


q* 


(3 = 6.0 


(3 = 6.2 


a = 


B 3 7 /2 


Unsmeared 
Unsmeared 


l/a 

n/a 


-1.951(09) 
0.606(06) 


-1.097(15) 
0.645(06) 




Smeared 
Smeared 


l/a 
n/a 


0.989(05) 
1.085(06) 


0.823(16) 
0.903(14) 


0.615(30) 
0.674(32) 


bT 


Unsmeared 
Unsmeared 


l/a 
n/a 


-1.486(11) 
0.822(05) 


-0.689(15) 
0.864(05) 




Smeared 
Smeared 


l/a 
n/a 


1.240(06) 
1.288(06) 


1.030(16) 
1.076(17) 


0.766(37) 
0.810(39) 



Table 3: Results for B? (NDR, 2 GeV) and Bf (NDR, 2 GeV), at the physical kaon mass, 
using a scales. 



contamination by excited states or off-shell matrix elements. Thus we do not give details 
concerning these other uncertainties. For our final results, we quote 

£? /2 (NDR, 2 GeV) = 0.62 ± 0.03(stat) ± 0.06(syst) , (34) 
£| /2 (NDR, 2 GeV) = 0.77 ± 0.04(stat) ± 0.04(syst) , (35) 

where the systematic error is our estimate (based on the q* variation) of the uncertainty due 
to the truncation of perturbation theory. 

It is interesting to compare our staggered results with those from Wilson fermions. While 
a continuum extrapolation is not yet available, the results at (3 = 6.0 || are 

Bj /2 (NDR, 2 GeV, Wilson) = 0.58 ± 0.02 ± 0.07 , (36) 
Ss /2 (NDR, 2 GeV, Wilson) = 0.81 ± 0.03 ± 0.03 . (37) 

The central values and error bars have been determined in the same way as in this paper. 
Obviously the agreement is already rather good, and indicates indirectly that the 0(a) errors 
in the Wilson case are not particularly large. This is in contrast to the case of Bk, where 
the explicit breaking of chiral symmetry disrupts the delicate cancellation between the VV 
and AA matrix elements, and gives dramatically large 0(a) errors. For B? and Bg the result 
is dominated by the PP matrix element, with no particular constraint from chiral symmetry, 
presumably yielding more moderate 0(a) errors. 
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6 Non-perturbative results for matching constants 



It is clear from the results of the previous two sections that using finite order perturbation 
theory to match lattice and continuum operators is an important source of uncertainty in 
results for matrix elements. In this section we investigate the accuracy of perturbative 
matching factors by calculating some of them non-perturbatively. In particular, we are able 
to assess the accuracy of the perturbative matching factor for the quark mass, Z m , which is 
an important ingredient in determining continuum quark masses. 

The basic idea is simple, and has been applied extensively with Wilson fermions. A given 
continuum operator can be discretized in different ways, each discretization having an asso- 
ciated matching factor. Only if these matching factors are chosen correctly will the different 
choices yield the same matrix elements. This allows a non-perturbative determination of the 
ratio of matching factors. Note that in such ratios the anomalous dimension factors cancel, 
implying that the ratios are finite functions of the lattice coupling. 

In this section, we apply this idea to the pseudoscalar and axial densities. It could, in 
principle, be applied also to four-fermion operators such as Qk, but this requires studying 
a matrix mixing problem, and thus using several external states, and is beyond the scope of 
the present work. We only consider operators for which the matching is diagonal. 

Our notation in this section differs from that used above. We use P UN , for example, to 
refer to the bare lattice operator constructed of unsmeared fields, 

>UN _ ft 1 /77~ u-untt q t \,,UN 



P UN = (l/sjN f )xr(K®ZM N - (38) 
In other words, we do not include the matching factor in the definition of P UN . The matching 



equation becomes [cf. Eq. fllq) 1 



pcont = Z UNpUN[ X + G ( a 2)] (39) 

Z UN = 1 + O^T) ( _ 7 (0) ln(g * a) + c UN } + 0(a 2 } (4Q) 
47T 

Here we have used the fact that for the matrix elements of interest the corrections are 
quadratic in the lattice spacing. Similar definitions apply to P SM ; A™ and A^ M . We also 
need to introduce the gauge- invariant version of the axial current A 1 ^ 1 , which is the unsmeared 
current with appropriate gauge links included. Note that P GI = P UN , i.e. the unsmeared 
pseudoscalar density is already gauge invariant since it is local. 

The quantities we determine non-perturbatively are Zj, M /Zp N , Z^ M /Z^ N , and . The 
first we obtain starting from the result 

(0\P cont \K°) = (0|Z™P UN |j^)[l + O(a 2 )] = (0|Zf M P SM |i^)[l + O(a 2 )] (41) 

from which we find the non-perturbative estimate 

Z SM / |pUN|jW)\ 

The equations for Z^/Za are identical except that P is replaced by A4. The one-loop 
perturbative results for the relevant ratios are 



Z 



SM 
P 



Z™ 



a- 



Ms(<f) 



1+ mow ' (d 
47T 



SM TJN\ 



a- 



1 + MbV * ' x 30.2532 , (43) 

47T 
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ySM 
vUN 



1 + 



a- 



Ms(?*) 



47T 



(4 X 



SM 



a- 



MSW* 



47T 



x 1.9384 



(44) 



The determination of 



proceeds slightly differently. We note that the gauge-invariant 
axial current is partially conserved on the lattice, implying that Z^ 1 = 1. Thus we could 
determine Z^( N by taking ratios of matrix elements of A Gl to those of While we have 

not calculated matrix elements using we do have available the matrix elements of its 
divergence, d^A Gl = 2mP Gl = 2mP . This is sufficient as long as the matrix element 
involves non-zero momentum transfer. Note that the lattice partial conservation equation is 
exact as long as we use the appropriate lattice derivative |23] 
arrive at 

1 f-sinh(m K )(0\A^ N \K^)^ 



Putting this all together we 



7UN 



2(m q /u )(Q\P™\K° G ) 



0(a z 



(45) 



where m q /uo is tadpole improved quark mass, and mx is the mass of the K G . All quantities 
on the r.h.s. of this equation are in lattice units. We have used sinh(m^) on the r.h.s., rather 
than rriK itself, because if we replace v4 UN by A Gl then the ratio on the r.h.s. is exactly equal 
to unity. In other words, sinh(m^) is the appropriate kinematical factor for the exactly 
conserved current. We choose to keep it for the unsmeared current in the hope that it will 
reduce the size of the 0(a 2 ) terms. The perturbative expression to which Eq. ( [45] ) should be 
compared is 



1 



1 + 



a- 



MS 



4:71 



7T~ 



9.17479) 



(46) 



Note that if we had used the average link in Landau gauge to determine uq, rather than the 
average plaquette, then Z™ = 1 at one-loop order. 

We determine the required ratios of matrix elements using the quantities previously used 
to determine the vacuum saturation approximants appearing in the S-parameters. For 
example, consider the ratio 



Rp(t) 



(Wjh) EyP VN (y, t))(Ey> P VN (tf, t)W(t 2 )) 
(W(ti)E^ SM (^))(E^ SM (^i)^(t 2 )) 



(47) 



For ti <gC t <§C t 2 this should be independent of t and gives directly (Zp M / Zp N ) 2 aside from 
0(a 2 ) corrections. We average Rp(t) over the same plateau regions as for the S-parameters. 
Similar ratios are used for the other quantities. 

The resulting data for Zp^/Zp N and Z^ M /Z^ N are well represented by a linear function 
of m 2 K . This dependence on m 2 K is an 0(p 2 a 2 ) discretization error, because any physical 
dependence cancels in the ratio of matrix elements, We remove this error by extrapolating 
to the chiral limit. We do a similar extrapolation for 1/Z^ N , although the dependence on 
m 2 K is much weaker, presumably because of the sinh(mx) factor in Eq. (53). 

Our non-perturbative results for the ratio of smeared to unsmeared matching factors 
after chiral extrapolation are collected in Table f|. What is most striking is the substantial 
dependence on lattice spacing, particularly for the ratio of Z P 's. This is due to a combination 
of 0(a 2 ) discretization errors and the variation of the perturbative matching factors which 
depend on g 2 (a). To analyze these results we assume the following form for the ratios 



RATIO(non-pert) = RATIO (one -loop; q*) + a 2 A 2 



(48) 
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Quantity 


Method 




(3 = 6.0 


(3 = 6.2 


a = 




Non-pert 




1.131(5) 


1.060(1) 


0.972(6) 


/Z^ N 


Pert(g* = 


l/a) 


0.970 


0.973 


0.977 




Pert(g* = 


it /a) 


0.979 


0.981 


0.983 




Non-pert 




2.245(11) 


1.861(6) 


1.380(19) 


Z|, M /Zp N 


Pert(g* = 


l/a) 


1.462 


1.416 


1.359 




Pert(g* = 


Tt/a) 


1.323 


1.301 


1.272 



Table 4: Non-perturbative and perturbative results for ratios of matching constants. 
Quadratic extrapolations to a = use a scales. 



Method 


(3 = 6.0 


(3 = 6.2 


(3 = 6.4 


Non-perturbative 


0.978(09) 


1.010(09) 


0.998(20) 


Pert (q* = 1 /a) 


0.986 


0.987 


0.988 


Pert (q* = vr/a) 


0.990 


0.991 


0.991 



Table 5: Results for 1/Z\ 



where the one-loop results are given above, and A is an unknown constant. In other words we 
ignore completely higher powers of a, and assume that higher powers of a are well represented 
by the appropriate choice of q* = K/a. The difference RATIO (non-pert)— RATIO (one- loop) 
should then, for the right choice of q*, extrapolate to zero in the continuum limit. Conversely, 
one could regard this procedure as providing an approximate non-perturbative definition of 
q*. To show the individual variations in the perturbative^] and non-perturbative results we 
give, in Table |], the continuum value for each obtained by linear extrapolation in a 2 . We do 
this for our two standard choices q*a = 1 and n. What we find remarkable is that, modulo 
the simplifying assumption of Eq. (f^8|), the non-perturbative and perturbative predictions 
agree if we use q* ~ l/a but not for q* pa it /a. The discrepancy for q* = n/a is particularly 
significant for Zf^/Zf*. 

The results for 1/Z^ N , given in Table [|, behave very differently. There is very little 
dependence on the lattice spacing — presumably because the unsmeared current is very sim- 
ilar to the gauge-invariant current for which all the numbers in the table would be unity 
independent of lattice spacing. In fact, the errors are such that an extrapolation to a = 
is not useful, and so we compare our results to perturbation theory at each lattice spacing. 
The results are reasonably consistent, but we cannot distinguish between different values of 
q* in this case. 

An important application of the above results is to estimate the reliability of the one-loop 

10 Note that we are here making use of the fact that our extrapolation to the continuum limit does not 
remove terms which vary logarithmically. This is an example of the problem discussed in sec. |j. 
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result for the matching factor Z m . This factor converts the lattice results for quark masses 
to a continuum scheme like MS, as discussed in Ref. The perturbative result, after 

tadpole improvement is 



= O = zW^t) = 1 " ^f^(-^ 0) ln (« + C ^ N ) > ( 49 ) 

where for simplicity we have chosen to consider the case where the final scale \x equals the 
matching scale q*. The one-loop correction to Z m is large at typical lattice spacings. For 
example, at (3 = 6, and taking q* = 1/a, Z m = 1.598. This suggests that higher order 
corrections may be important. We can, however, rewrite Z m as 

Zm= \z^)\zW) ■ (50) 

The results of Table |] show that the bulk of the perturbative correction lies in the first factor, 
Zf M /Z£ N . At p = 6 it is 1.462, while the second factor is 1.136 (again for p = q* = 1/a). 
Thus, one would expect that the dominant source of higher order terms in Z m is the first 
factor, and that the uncertainty which they introduce could be substantially reduced by 
obtaining a non-perturbative estimate of this factor. We have attempted such an estimate 
above, with the preliminary conclusion that perturbation theory with q* = 1/a works to 
within a few percent, aside from discretization errors. If we accept this result then we obtain 
the partly non-perturbative estimate^ Z m = 1.462 x 1.136 = 1.66. The point we wish to 
stress is that this is not very different from the one-loop estimate of 1.60. In particular, the 
difference is much smaller than the naive estimate of the two-loop contribution, 0.6 2 = 0.36, 
based on the assumption of geometric growth. This analysis thus suggests that the one-loop 
perturbative value of Z m used in the analysis of quark masses is good to about 5% for (3 > 6.0 
provided one uses o>y^(q* ~ 1/a). 



7 Conclusions 

In this paper we have presented a variety of results for weak matrix elements using staggered 
fermions. Our major focus has been on the importance of using a variety of discretizations 
of continuum operators. There are two reasons for doing so. First, comparing results with 
different lattice operators gives an estimate of the uncertainty in the matching factors be- 
tween continuum and lattice operators. For this may be the dominant source of error 
in future calculations, aside from that due to quenching. Second, for some operators the 
perturbative matching factors are not convergent at present couplings, and so one must use 
different discretizations. It turns out that the smeared operators have uniformly moderate 
perturbative corrections. Using them we are able to obtain the first results for B^ 2 and B^ 2 
using staggered fermions. 

Our results for the 5-parameters confirm and extend existing lattice results. In particular, 
for Br we find that smeared operators give results consistent with those from unsmeared 
and gauge-invariant operators. We confirm the low value found in our preliminary study 0, 

11 It is not advantageous to directly use the non-perturbative results, e.g. Zf, M /Zp N = 2.245 at f3 = 6. 
Doing so introduces additional 0(a 2 ) errors which one would then have to remove by extrapolation. 
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a result which has been improved and extended by the JLQCD collaboration 0. For B^ 2 

3 /2 

and B 8 we find results consistent with those using Wilson fermions. All these numbers are 
important inputs into analyses attempting to constrain the CKM matrix. It is encouraging 
that the errors we are considering are at the few percent level. It is important to stress, 
however, that we are still using the quenched approximation, and also working with a kaon 
composed of degenerate quarks. For a discussion of the importance of these approximations 
see Refs. [|T], [13]. 



As an offshoot of our study, we have calculated several ratios of matching factors non- 
perturbatively. These ratios are finite functions of the lattice coupling, and thus allow a 
test of tadpole improved perturbation theory. We find that one-loop perturbation theory 
works well if we set the scale in the one-loop coupling, a^fg(g*), to be q* = 1/a, but not 
for q* = 7c/a. This is consistent with the expectations of Ref. ||16|| . This conclusion is, 
however, preliminary because we have results only at two lattice spacings. To convincingly 
disentangle discretization errors from perturbative corrections will require precise results at 
several lattice spacings. 

We have used the results for ratios of matching factors to make a partly non-perturbative 
estimate of the size of the matching factor for the quark mass, Z m . Our result is ~ 5% higher 
than the one- loop perturbative result. This is a small enough change that it does not alter 



the essential conclusion of Ref. ||12|| , namely that light quark masses are considerably smaller 
than previously thought. 

Finally, we note that our results show many examples of significant discretization errors. 
To make progress with simulations of full QCD, where one is restricted to larger lattice 
spacings, it may be necessary to improve the staggered fermion action. 
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